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DD Correlations as a Sensitive Probe for Thermalization in High-Energy Nuclear Collisions 
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We propose to measure azimuthal correlations of heavy-flavor hadrons to address the status of thermalization 
at the partonic stage of light quarks and gluons in high-energy nuclear collisions. In particular, we show that 
hadronic interactions at the late stage cannot significantly disturb the initial back-to-back azimuthal correlations 
of DD pairs. Thus, a decrease or the complete absence of these initial correlations does indicate frequent 
interactions of heavy-flavor quarks and also light partons in the partonic stage, which are essential for the early 
thermalization of light partons. 
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I. INTRODUCTION 

Lattice QCD calculations, at vanishing or finite net-baryon 
density, predict a cross-over transition from the decon- 
fined thermalized partonic matter (the Quark Gluon Plasma, 
QGP) to hadronic matter at a critical temperature T c w 150— 
180MeVd. 

Measurements of hadron yields in the intermediate and high 
transverse momentum (pj) region indicate that dense matter 
is produced in Au+Au collisions at RHIC |2|, |3[. The ex- 
perimentally observed large amount of elliptic flow of multi- 
strange hadrons |4], such as the (j) meson and the Q. baryon, 
suggest that collective motion develops in the early partonic 
stage of the matter produced in these collisions. A crucial is- 
sue to be addressed next is the thermalization status of this 
partonic matter. 

Heavy-flavor (c, b) quarks are particularly excellent 
tools ||5l la, 0] to study the thermalization of the initially cre- 
ated matter. As shown in Fig.^ their large masses are almost 
exclusively generated through their coupling to the Higgs field 
in the electro-weak sector, while masses of light quarks (u, 
d, s) are dominated by spontaneous breaking of chiral sym- 
metry in QCD. This means that in a QGP, where chiral sym- 
metry might be restored, light quarks are left with their bare 
current masses while heavy-flavor quarks remain heavy. Due 
to their large masses (3> AgcoX the heavy quarks can only 
be pair-created in early stage pQCD processes. Furthermore, 
their production cross sections in nuclear collisions are found 
to scale with the number of binary nucleon-nucleon colli- 
sions iH 0. In the subsequent evolution of the medium, 
the number of heavy quarks is conserved because the typical 
temperature of the medium is much smaller than the heavy 
quark (c, b) masses, resulting in negligible secondary pair 
production. In addition, the heavy quarks live much longer 
than the lifetime of the formed high-density medium, de- 
caying well outside. These heavy quarks (c, b) can partic- 
ipate in collective motion or even kinetic ally equilibrate if, 
and only if, interactions at the partonic level occur at high 
frequency. The idea of statistical hadronization of kineti- 
cally equilibrated charm quarks Ml Oil even predicted signifi- 



cant changes in hidden charm hadron production 1 1 lj - Hence, 
heavy-flavor quarks are an ideal probe to study early dynamics 
in high-energy nuclear collisions. Recent STAR and PHENIX 
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FIG. 1: (Color online) Quark masses in the QCD vacuum and the 
Higgs vacuum. A large fraction of the light quark masses is due to 
the chiral symmetry breaking in the QCD vacuum. The numerical 
values were taken from Ref. 1 1 2ll . 



results on elliptic flow and nuclear modification factors of 
non-photonic electrons indicate that charm quarks might in- 
deed participate in the collective motion of the matter pro- 
duced in Au+Au collisions at RHIC. To explain the data, a 
large drag diffusion coefficient in the Langevin equation, de- 
scribing the propagation of charm quarks in the medium, is 
found to be necessary fl3ll . Two- and three-body interac- 
tions lfl4l ll 511 of heavy-quarks and resonant rescattering itldl 
in the partonic stage seem to become important. These inves- 
tigations suggest that heavy-quarks actively participate in the 
partonic stage. 

In this paper, we study the change of azimuthal correla- 
tions of D and D meson pairs as a sensitive indicator of fre- 
quent occurrences of partonic scattering. Since heavy-flavor 
quarks are pair-created by initial scattering processes, such as 
gg — > cc, each quark-antiquarkpair is correlated in relative az- 
imuth A(j) due to momentum conservation. In elementary col- 
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lisions, these correlations survive the fragmentation process 
to a large extent and hence are observable in the distribution 
of the relative azimuth of pairs of D and D mesons fnll . We 
evaluate by how much these correlations should be affected by 
the early QGP stage and by the hadronic scattering processes 
in the late hadronic stage. 



II. RESULTS AND DISCUSSIONS 

To explore how the QCD medium generated in central ultra- 
relativistic nucleus-nucleus collisions influences the correla- 
tions of D and D meson pairs, we employ a non-relativistic 
Langevin approach which describes the random walk of 
charm quarks in a QGP as was first described in Refs. (USE! . 
Here, both the drag coefficient y and the momentum-space dif- 
fusion coefficient a depend on the local temperature, T with 
a momentum-independent y(r) = aT 2 . We chose a — 2 ■ 10~ 6 
(fm/c)- 1 MeV" 2 . Details of our calculations can be found 
in CI. 

In order to isolate the effects purely due to parton-parton 
re-scattering in the outlined medium, we generated c and 
c quarks with the same pj and zero longitudinal momen- 
tum back-to-back, i.e. with A(j) = n at a time of the order 
of 1 / m c ~ 0. 1 fm/c, and used a delta function for fragment- 
ing the charm quark into a charmed hadron at the hadroniza- 
tion stage. The radius r where each pair is created is ran- 
domly generated from a distribution reflecting the number 
of binary nucleon-nucleon collisions taking place at that ra- 
dius, p(r)dr °= (R 2 — r 2 )2nrdr. The angular distribution of 
the initially produced charm quarks in the transverse plane is 
isotropic. 

The evolution of the charm momenta with the Langevin 
equation is stopped when T reaches the critical temperature 
T c = 165 MeV or when the charm quark leaves the QGP vol- 
ume. Furthermore, to get a first estimation of the QGP ef- 
fect on the charm quark pairs azimuthal correlation, we omit 
the possible contribution of the mixed phase. Figure 0(a) 
shows the results for the D meson (charm quark) pairs an- 
gular correlations for different initial charm quark pj values, 
for 7b = 300 MeV and x = 0.5 fm/c (typical values for RHIC 
collisions). We see that the fastest charm quarks (represented 
by the pj — 3 GeV/c) are able to escape from the QGP without 
suffering significant medium effects, while the slower quarks 
(see the pj — 0.5 GeV/c black line) have their pair azimuthal 
correlation almost completely smeared out by the interactions 
in the medium. 

Figure 0(b) shows the corresponding results for To = 
700 MeV and To = 0.2 fm/c, values representative of LHC 
energies. Here, the interactions of the charm quarks with the 
medium are so frequent that only the most energetic charm 
quarks preserve part of their initial angular correlation; low pj 
pairs can even be completely stopped by the medium. Because 
the c and the c quarks of a given pair are created together, in 
the same space point, the pair has a higher escape probability 
if both quarks escape the medium from the side of the fireball 
where the thickness is smaller. Thus, the A(j) distribution is 
shifted towards smaller values. 



Presently, reliable numbers on the drag coefficient do not 
exist lfl6l fl9ll . The dependence on the drag coefficient is 
shown in Figs. |2](c) and (d) for the cc angular correlations 
for energetic c quarks: pj = 3 and 10 GeV/c for To = 300 and 
700 MeV, respectively. These correlations vanish when a is 
increased by around a factor of five in the first case and by 
around a factor of two in the latter one, with respect to values 
from pQCD-calculations. 
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FIG. 2: (Color online) Correlations in relative azimuth A(|) of DD 
pairs from Langevin calculations with Tq = 300 MeV (RHIC) and 
700 MeV (LHC). The upper part shows the dependence of the corre- 
lations on the initial pj of the c quark; the lower part shows the drag 
coefficient dependence. 

We have demonstrated that initial correlations of cc pair 
correlations are clearly affected by the formed QCD medium, 
and sensitive to its temperature and drag coefficient. The 
question is to what extent hadronic interaction can modify 
these initial correlations. In the following, we apply the mi- 
croscopic hadronic transport approach UrQMD v. 2.2 |'2(X EH. 
A phase transition to a QGP state is not incorporated into the 
model dynamics. In the latest version of the model, PYTHIA 
(v. 6.139) was integrated to describe the energetic primary el- 
ementary collisions 1 22]. Measured particle yields and spectra 
are well reproduced by this approach 12211. 

In this model, D mesons stem from the early-stage high- 
energy nucleon-nucleon collisions, calculated with PYTHIA. 
For their propagation in the hadronic medium, we consider 
elastic scattering of D mesons with all other hadrons. The 
hadronic scattering cross-section for D mesons is generally 
considered to be small 1 23] and we take 2 mb in our cal- 
culation. The results of the UrQMD calculations are shown 
in Fig. |3](a), for minimum bias Au+Au reactions at ^/s NN = 
200 GeV To better illustrate the change in the angular corre- 
lations of the DD pairs, we show the ratio between the final 
distribution, affected by the evolution with UrQMD, and the 
initial one, directly obtained from PYTHIA. In case hadronic 
interactions do not modify the correlations, this ratio is unity. 

For this specific analysis, all DD pairs were selected, ir- 
respective of the pj of the D mesons. Fig. |3](a) shows our 
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results assuming a forward peaked scattering cross section, 
exp(7 cos 0), where is the scattering angle in the center of 
mass system of the decaying resonance. Even for a large cross 
section of 20 mb, the DD pair correlations are barely affected 
by hadronic interactions. 

On the other hand, intermediate resonances formation (es- 
sentially D+7I <->D*), might result in isotropic emission, i.e. 
the correlation might be completely destroyed by a single scat- 
tering process. Only when assuming the unlikely case of a 
large cross section of 20 mb with fully isotropic emission, the 
DD pair correlations are significantly modified, see Fig. 0(b). 
We also find that most scatterings occur at relatively early 
times, when the hadron density is high. At RHIC energies, 
the (pure) hadronic stage is presumably shorter than assumed 
in the present calculations. 

Hence, we conclude that hadronic interactions are un- 
likely to significantly modify DD pair correlations. A change 
in these angular correlation must be dominated by frequent 
parton-parton scatterings occurring in the QGP phase. 
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FIG. 3: Ratios between the A<|> distributions of DD pairs produced in 
minimum bias Au+Au reactions at ^/s NN = 200 GeV, before and after 
hadronic rescattering, using forward (left panel) and isotropic (right 
panel) angular distributions. 

The calculations reported above show that heavy-flavor cor- 
relations provide a sensitive tool to directly access the sta- 
tus of equilibration of the partonic medium. A similar pic- 



ture emerges from the observation and suppression of back-to- 
back correlations of unidentified high pj hadrons (jets). Both 
features originate from the propagation of partons (heavy 
quarks or jets) in the medium. However, D mesons (and B 
mesons) have the advantage that they can still be identified, 
even if they lose a significant fraction of their energy and get 
kinetically equilibrated. 



III. CONCLUSION AND OUTLOOK 

In summary, we argue that the observation of broadened an- 
gular correlations of heavy-flavor hadron pairs in high-energy 
heavy-ion collisions would be an indication of thermalization 
at the partonic stage (among light quarks and gluons). We 
have seen that hadronic interactions at a late stage in the colli- 
sion evolution cannot significantly disturb the azimuthal cor- 
relations of DD pairs. Thus, a visible decrease or the complete 
absence of such correlations, would indicate frequent interac- 
tions of heavy-flavor quarks and other light partons in the par- 
tonic stage, implying early thermalization of light quarks in 
nucleus-nucleus collisions at RHIC and LHC. 

These measurements require good statistics of events in 
which both D mesons are cleanly reconstructed. A complete 
reconstruction of the D mesons (i.e. of all their decay prod- 
ucts) in full azimuth is essential to preserve the kinematic in- 
formation and to optimize the acceptance for detecting corre- 
lated D meson pairs. Solid experimental measurements in pp 
and light-ion collisions, at the same energy, are crucial for de- 
tailed studies of changes in these azimuthal correlations, and 
should be per formed as a function of pj. Proposed upgrades 
of STAR (3 and PHENIX at RHIC and the ALICE detec- 
tor ll25ll at LHC with micro-vertex capabilities and direct open 
charm reconstruction should make these measurements possi- 
ble. 
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